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Abstract 

We derive new constraints on the quantities Syr Y i X,Y = L,R, which parametrise the flavour-off- 
diagonal terms of the charged slepton mass matrix in the MSSM. Considering mass and anomalous 
magnetic moment of the electron we obtain the bound I^^^r^I ^ 0. 1 for tan/3 = 50, which involves 
the poorly constrained element S]^ R . We improve the predictions for the decays r — > fij, r — > 
e7 and /i — > by including two-loop corrections which are enhanced if tan/3 is large. The finite 
renormalisation of the PMNS matrix from soft SUSY-breaking terms is derived and applied to the 
charged-Higgs-lepton vertex. We find that the experimental bound on BR(t — > ej) severely limits 
the size of the MSSM loop correction to the PMNS element U e 3, which is important for the proper 
interpretation of a future U e ^ measurement. Subsequently we confront our new values for 8 l [ L with a 
GUT analysis. Further, we include the effects of dimension-5 Yukawa terms, which are needed to fix 
the Yukawa unification of the first two generations. If universal supersymmetry breaking occurs above 
the GUT scale, we find the flavour structure of the dimension-5 Yukawa couplings tightly constrained 
by (J. -t ej. 

1. Introduction 

Weak-scale supersymmetry (SUSY) is an attractive framework for physics beyond the standard model (SM) of 
particle physics. The SM fields are promoted to superfields, with additional constituents of opposite spin. Due 
to their identical couplings, they cancel the quadratic divergent corrections to the Higgs mass. Since none of the 
SUSY partners have been observed in experiments, supersymmetry must be broken and the masses of the SUSY 
partners are expected to be in the multi-GeV region. 

A supersymmetric version of the standard model requires a second Higgs doublet in order to cancel the Higgsino- 
related anomalies and to achieve electroweak symmetry breaking. At tree level, one of the Higgs doublets, H u , 
couples to the up-type particles, whereas the other doublet, Hj, couples to the down type particles. The Yukawa 
couplings of the minimal supersymmetric standard model (MSSM) read 

IUmssm = Y* j u1QjH u + Y^dtQjHa + >/'V/.,//./ + ^H d H u . (la) 
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Neutrinos are massless in the MSSM; however, experiments and cosmological observations consistently point at 
small but non-vanishing masses in the sub-eV region. We will therefore consider an extended MSSM with three 
right handed neutrinos, where the Yukawa couplings are given by 

W = Wmssm + ufLjH u + ^Mp?v] . (lb) 

Here, Q and L denote the chiral superfields of the quark and lepton doublets and u c , d c , e c and v c the up and 
down-quark, electron and neutrino singlets, respectively. Each chiral superfield consists of a fermion and its scalar 
partner, the sfermion. The Yukawa coupling matrices Y u ,d,i,v are defined with the right-left (RL) convention. The 
field v c is sterile under the SM group, so we allow for a Majorana mass term in addition to the Dirac coupling. 
The respective mass matrix is denoted by Mr and the scale of Mr is undetermined but expected to be above the 
electroweak scale, M ew (see Sec. [J). 

The Higgs fields acquire the vacuum expectation values (vevs) 

(H u ) = v u , (H d ) = v d . (2) 

where \v u \ 2 + \v d \ 2 = v 2 = (174 GeV) 2 . The ratio of the two vevs is undetermined and defines the parameter 
tan/3, 

— =:tan/3. (3) 

Vd 

While tan j3 is a free parameter of the theory, there exist lower and upper bounds on its value. Experimentally, 
Higgs searches at LEP rule out the low-tan /3 region in simple SUSY models |Tj . This result fits nicely with the 
theoretical expectation that the top Yukawa coupling should not be larger than one. The region of the MSSM 
parameter space with large values of tan /3 is of special importance for the flavour physics of quarks and leptons. 
We therefore have a brief critical look at the upper bounds on this parameter: Demanding a perturbative bottom 
Yukawa coupling y^ naively leads to an upper limit on tan /3 of about 50 inferred from the tree-level relation 

m b m b 

Vb = a ~ tan/3 . 4 

v ■ cos p v 

Similarly, the MSSM provides a natural radiative breaking mechanism of the electroweak symmetry as long as 
i/b < yt at a low scale [2] ■ At tree level, the ratio of the Yukawa couplings is given by 



1 > 



^tan/3. (5) 
m t 



Since m t ([i)/mb(fi) ~ 60 at the electroweak scale, tan/3 should not exceed this value. 

Both arguments, however, do only hold at tree level. In particular, down quarks as well as charged leptons 
couple to H u via loops. As a result, if we take tan /3-enhanced contributions into account, an explicit mass 
renormalisation changes the relation of Yukawa coupling and mass [SHE] . The tan /3 enhancement of the bottom 
coupling in Eq. ^ can be compensated; similarly, the ratio of the Yukawa couplings is changed due to an explicit 
bottom quark mass renormalisation. We will find that values of tan/3 up to 100 both provide small enough 
Yukawa couplings and do not destroy natural electroweak symmetry breaking. 

Large values for tan/3 are interesting for two reasons. One, in various grand-unified theories (GUTs), top and 
bottom Yukawa couplings are unified at a high scale. In this case, it is natural to expect tan /3 = mt/rrib, as shown 
above. Two, many supersymmetric loop processes are tan /3-enhanced due to chirality-flipping loop processes with 
supersymmetric particles in the loop. This enhancement can compensate the loop suppression and therefore large 
values of tan/3 lead to significant SUSY corrections. 

In this paper, we will study the lepton sector in the (extended) MSSM. Since the neutrinos are massive, the 
leptonic mixing matrix, t/pMNS? is n0 longer trivial and leads to lepton flavour violation (LFV). In its standard 
parametrisation, it reads 




C/PMNS = C 23 S23 1 s 12 c 12 e 4 ^ , (6) 
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with Sij = sin (9^ and cy = co&Oij. The two phases cti^ appear if neutrinos are Majorana particles. They are 
only measurable in processes which uncover the Majorana nature of neutrinos, such as neutrinoless double beta 
decay. 

The PMNS matrix allows for flavour transitions in the lepton sector, in particular neutrino oscillations, through 
which its parameters are well constrained. Compared with the mixing angles of the quark mixing matrix, Vckm, 
two mixing angles, namely the atmospheric and solar mixing angles, # 2 3 = $atm an d #12 = 9 m \, are surprisingly 
large, whereas the third mixing angle is small. The current experimental status at la level is as follows 00 

9 l2 = 34. 5 ± 1.4 , Amlx = 7. 67±°;^ ■ 1(T 5 eV 2 , 

-2. 37 ± 0. 15 • 1CT 3 eV 2 inverted hierarchy, 
2. 46 ± 0. 15 • 10 -3 eV 2 normal hierarchy, 

? l3 = 0.0t™. (7) 



7 23 = 42.3lii, Amtx = <i io Ac ^ n ic in _ 3 Tr2 



These values are determined by the atmospheric and solar mass splitting Am 2 tm = Am 2 3 , Am 2 ol = Am^, leaving 
the absolute mass scale open. The pattern of mixing angles is close to tri-bimaximal, corresponding to 623 = 45°, 
9x2 — 35°, and 9x3 = 0° [8j. Due to the smallness of 6x3, the CP phase S is unconstrained. Tri-bimaximal mixing 
can be motivated by symmetries (see Ref. [9] and references therein), which constrain fundamental quantities 
like Yukawa couplings or soft SUSY-breaking terms. Measurable quantities like U e 3 usually do not point directly 
to fundamental parameters, but are sensitive to corrections from all sectors of the theory. The analysis of such 
corrections is therefore worthwhile. A large portion of this paper is devoted to the influence of supersymmetric 
loops and higher-dimensional Yukawa terms on observables in the lepton sector of the MSSM. 

In a supersymmetric framework, additional lepton flavour violation can be induced by off-diagonal entries in the 
slepton mass matrix, which parametrise the lepton-slepton misalignment in a model independent way. However, a 
generic structure of the soft masses is already excluded because too large decay rates for lj — > ^7 would arise. To 
avoid this flavour problem, the SUSY breaking mechanism is often assumed to be flavour blind, yielding universal 
soft masses at a high scale. Then the PMNS matrix is the only source of flavour violation in the lepton sector, as 
is the CKM matrix for the quarks; this ansatz is called minimal flavour violation. The soft terms do not cause 
additional flavour violation and the various mass and coupling matrices are flavour-diagonal at some scale in the 
basis of fermion mass eigenstates, e.g., 



- m\ = ml 1 , m\ = m\ = , M = A Yi . (8) 



L 



Here, m 2 - . denote the soft mass matrices of the sleptons (see Eq. (|87p). m 2 H the analogous soft masses of the 
Higgs doublets, and Ai is the trilinear coupling matrix of the leptons. 

Even if the soft terms are universal at the high scale, renormalisation group equations (RGE) can induce non 
vanishing off-diagonal entries in the slepton mass matrix at the electroweak scale. Lepton flavour violation can be 
parametrised by non-vanishing S l ^ Y a t t ne electroweak scale in a model-independent way, where 5 l -^ Y is defined 
as the ratio of the flavor-violating elements of the slepton mass matrix (f8"7| and an average slepton mass (see 
Eq. Ull)), 

S£ Y = ; gl= , X,Y = L,R, i,j = 1,2,3 (i^j). (9) 

m 2 ix m 2 jY 

The flavour-off-diagonal elements Am l ^- Y are defined in a weak basis in which the lepton Yukawa matrix Yi in 
Eq. {Taj is diagonal. According to the chiralities of the sfermion involved, there are four different types, 8ll, Srr, 



8lr, and Srl- The tolerated deviation from alignment can be quantified by upper bounds on S^y, as discussed 
above and are already extensively studied in the literature (see for example [lUffTS] and references therein). 

Being generically small, the sfermion propagator can be expanded in terms of these off-diagonal elements, 
corresponding to the mass insertion approximation (MIA) |13[fT4] . Instead of diagonalising the full slepton mass 
matrix and dealing with mass eigenstates and rotation matrices at the vertices, in MIA one faces flavour-diagonal 
couplings and LFV appears as a mass insertion in the slepton propagator. This approach is valid as long as 
^xy ^ 1 an d makes it possible to identify certain contributions easily. For a numerical analysis an exact 



1 Recently, a hint for non-zero 613, sin 2 $13 = 0. 016 ± 0. 010 (lcr), was claimed in Ref. [7|. 
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diagonalisation of all mass matrices is, of course, possible. In Ref. |12) a systematic comparison between the full 
computation and the MIA both in the slepton and chargino/neutralino sector clarifies the applicability of these 
approximations. 

This paper provides a comprehensive analysis of the lepton sector in the MSSM, focusing on the phenomeno- 
logical constraints on the parameters S l ^ Y m Eq. ©■ In Sec. 2 we briefly review the supersymmetric threshold 
corrections to in Eq. ((4]) and relax the usually quoted upper bounds on tan j3. In Sec. 3 we derive new con- 
straints on the <^y' s by studying loop corrections to the electron mass, finite renormalisations of the PMNS 
matrix and the magnetic moment of the electron at large tan /3. As a byproduct we identify all tan /3-enhanced 
corrections to the charged-Higgs coupling to leptons. We then improve the MSSM predictions for the decay rates 
of lj — > ij7 by including tan /3-enhanced two-loop corrections. In Sec. 4 we embed the MSSM into SO(10) GUT 
scenarios and study RGE effects. Even for flavour-universal soft breaking terms at the GUT scale Mqut sizable 
flavour violation can be generated between Mqut and the mass scales of the right-handed neutrinos |15| . Com- 
paring the results to the upper bounds on |<5^fI found in Sec. 3 enables us to draw general conclusions on GUT 
scenarios. We include corrections to the Yukawa couplings from dimension- 5 terms and constrain their possible 
flavour misalignment with the dimension-4 Yukawa matrices. In Sec. 5 we summarise our results. Our notations 
and conventions are listed in three appendices. 

2. Upper bound for tan (3 

The tree level mass of a particle receives corrections due to virtual processes. Supersymmetric loop corrections 
are small unless a large value for tan /3 compensates for the loop suppression. We will therefore start with a 
discussion of tan /3-enhanced loops and derive a relation between Yukawa coupling and mass, coming from the 
resummation of tan /3-enhanced corrections to the mass to all orders. 

Corrections to the mass with more than one loop and more than one coupling to Higgs fields do not produce 
further factors of tan /3. Nevertheless, tan /3-enhanced loops can become important in an explicit mass renor- 
malisation, where tan /3-enhanced contributions to all orders are taken into account, because counterterms are 
themselves tan /3-enhanced. 

Down-quarks and charged leptons receive tan /3-enhanced corrections to their masses, due to loops with H u . 
As a coupling to H u does not exist at tree level, tan /3-enhanced loops are finite; there is no counterterm to this 
loop induced coupling. Moreover, tan /3-enhanced contributions to self-energies do not decouple. The coupling of 
H u to the charged slepton is proportional to \x with /i = ©(Afgusy)- On the other hand, the integration over the 
loop momentum gives a factor 1/Msusy- Thus the dependence on the SUSY mass scale cancels out. In the large 
tan /3 regime, neutralino-slepton and chargino-sneutrino loops can significantly change the relation between the 
Yukawa couplings and masses [PB1 - I18) . 

oo phys 

m< 0) = mf ys + £ (-A)" mf hys = , (10) 

n— 1 

where the corrections A; are related to the self-energy £/ as A; = — S;/m;. This relation includes all tan/3- 
enhanced contributions and can be determined by only calculating two diagrams, according to chargino-sneutrino 
and neutralino-slepton loops (Figure [T]). As a result, the relation between Yukawa coupling and physical mass is 
given by 

(0) phys 

The individual contributions from the two diagrams to the self-energy S = S^* + T,*° are 

= 1^2 E X>> A I? 7 'ir' /'/ />V. (M| ?) m?) . (12) 

1=1,2 ] = l 
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Figure 1: Contribution to the self-energy E arising from chargino-sneutrino and neutralino-slepton loops 




Figure 2: Bottom quark coupling as a function of tan /? for eb — (black-dashed) and t\, — 0. 008 (red-solid). 

The loop integrals, couplings and rotation matrices are defined in the appendices lAl and iBl M^± and M«o denote 
the chargino and neutralino masses, respectively. The tan /3-enhanced transitions require a chirality flip and are 
given by 



A; = = ei tan /3 = — Mi [i tan/3 

mi 4-7T 



1 



-h \Mi^\m\\ - h (Mi,^,mi) + h [Mf,mf L ,mf R 



Q-2 

— M 2 li tan/3 

47T 



l -h (M|, /i 2 , m?J + f x (M|, Li 2 , ml) 



(13) 



with the loop- function fi(x,y,z) defined in Eq. (|74|) . 

The improved relation between the Yukawa coupling and the physical mass relaxes the upper bound for 
tan j3, as indicated in the Introduction. Equation (HJ changes to [TB] 



Vb 



mb tan f3 



v 1 + £b tan j3 

For down-quarks, the SUSY contributions are dominated by gluino loops such that 

£b - -5— mg Ll fx [T 
07T 



(14) 



(15) 



A typical value is et s» 0.008, leading to yj, = 0(1) for tan/3 ps 100 (see Fig. [2]). Similarly, the bound for natural 
electroweak symmetry breaking shifts to tan/3 < 100. 

In the following, we will use this relation to constrain lepton flavour violating parameters. 



3. Constraints on the flavour-violating parameters 

Various processes can be used to constrain lepton flavour violating (LFV) parameters. Remarkably, we can also 
use lepton flavor conserving (LFC) observables, due to double lepton flavour violation (LFV). Two LFV vertices 
lead to lepton flavor conservation (LFC) and so contribute to the LFC self-energies. In a similar manner, we will 
consider multiple flavour changes contributing to the magnetic moment of the electron. In addition, we consider 
LFV processes, in particular radiative decays. 
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As mentioned in the Introduction, we introduce the dimensionless parameters 8 l ^ Y YVA 

Am^ Y = 4 Y (rff^y) 2 = 6% Y y/mlml , X,Y = R,L, i^j. (16) 

Note that Amjjy has mass-dimension two. Both mf and Arn-^ Y are the diagonal and off-diagonal entries of 
the slepton mass matrix (|57|) . so rn^- Y is an average slepton mass. 

The effects discussed in this section stem from chirally enhanced self-energies, which involve an extra factor of 
tan/3 compared to the tree- level result, analogous to et in Eq. (fT5|) . Such effects have been widely studied before 
in the quark sector, yet most authors have performed their studies in the decoupling limit Msusy 3> v, where 
Msusy denotes the mass scale of the soft SUSY breaking parameters. (For a guide through the literature see 
Ref. [TS].) If one relaxes the condition Msusy 3> v, novel effects (namely those which vanish like some power of 
w/Msusy) can be analysed. Analytical results covering the case Msusy ~ v have been derived in Refs. |16lfl8H20] . 
numerical approaches were pursued in Refs. |21U22| . Superparticle contributions to physical processes vanish for 
Msusy — >• oo, typically as v 2 /M| USY , and can only be addressed with the methods of the latter papers. However, 
by combining the decoupling supersymmetric loop with resummation formulae valid for Msusy 3> v one can 
correctly reproduce the resummed all-order result to leading non- vanishing order in v/ Msusy- This approach has 
been used in an analysis of "flavoured" electric dipole moments in Ref. [23] . 

The possibility to constrain S l ^- Y through LFC processes has been pointed out in Ref. [24], which addresses 
leptonic Kaon decays. Here we analyse the constraints from two LFC observables which have not been considered 
before: In Sec. 13.11 we apply a naturalness argument to the electron mass, demanding that the supersymmetric 
loop corrections are smaller than the measured value. In Sec. 13.51 we study the anomalous magnetic moment of 
the electron. 



3.1. Flavour-conserving self-energies 

The masses of the SM fermions are protected from radiative corrections by the chiral symmetry "J — > e 1 " 75 ^. 
According to Weinberg, Susskind and 't Hooft, a theory with small parameters is natural if the symmetry is 
enhanced when these parameters vanish. The smallness of the parameters is then protected against large radiative 
corrections by the concerned symmetry. This naturalness principle makes the smallness of the electron mass 
natural. Radiative corrections are proportional to the electron mass itself Sm e oc m e ln(A/m e ) and vanish for 
m e = 0. If such a small parameter is composed of some different terms and one does not want any form of fine- 
tuning, one should require that all contributions should be roughly of the same order of magnitude; no accidental 
cancellation between different terms should occur. Hence, the counterterm of the electron mass should be less 
than the measured electron mass, 



dm e 




m e 0) - 


phvs 

rrie 




phvs 

rrie 



This naturalness argument for the light fermion masses was already discussed in the pioneering study of Ref. |10| . 
Since the authors of this analysis want to provide a model-independent analysis on classes of SUSY theories they 
consider only photinos and do not include flavour violation in the corrections to light fermion masses. Their 
derived upper bound for 5}} R depends on the overall SUSY mass scale and actually becomes stronger for larger 
SUSY masses. The case of radiatively generated fermion masses via soft trilinear terms is studied in [25] and 
an updated version including two flavour-violating self-energies can be found in |26| . Here we concentrate on 
the chirality-conserving flavour- violating mass insertions and use this argument to restrict the product 5j^ L S R 3 R 
which is then independent of the SUSY scale. Considering double lepton flavour violation, we demand that the 
radiative corrections should not exceed the tree-level contribution. For the electron mass, the dominant diagrams 
involve couplings to the third generation. As a result, we can constrain the product | ^^Jjyj . Note that \S RR \ 
has so far been unconstrained from radiative decay lj — > lij as we will discuss shortly in the following Section [3. 21 
This cancellation of the RR sensitivity is analysed in [TT]|T2] with the conclusion that even a better experimental 
sensitivity on BR(lj — > ^7) can not help to set strong constraints in the RR sector. However, with double mass 
insertions and the bound from \i — ¥ e^f it is possible to derive bounds on products like Sj^ L 6 RR . 

The diagram in Figure [3ja) achieves an m T tan (3 enhancement only if there is a helicity flip in the stau 
propagator. Since a\j (4tt) 3> Y" e 2 /(l67r 2 ), the higgsino contribution is negligible. A chargino loop can also 
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Figure 3: (a) LFC self-energy through double LFV and (b) dominant double LFV contribution to the electron 
mass renormalisation. 
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Table 1 : Snowmass Points and Slopes [27] f and two additional scenarios (masses in GeV) 



be neglected, because only left-left (LL) insertions for the sneutrinos can be performed and the helicity flip is 
associated with an electron Yukawa coupling. The left-right (LR) insertions are either not associated with an 
tan /3-enhanced contribution or suppressed by v/A/susy, compared to right-right (RR) and LL-insertions. Thus 
the dominant diagram involves a Bino and the scalar tau-Higgs coupling, as shown in Figure [3fb). For simplicity, 
we choose all parameters real and obtain 



For equal SUSY masses this simplifies to v = ~ jf^: m {l^ ^ $rr ■ This term is proportional to m T , in contrast 
to the LFC self energy, which is proportional to m e . Thus the counterterm receives an additional constant term 



— i?nf hys = — mf hys Ai *^f V —i6rni —imf^ 11 

Substituting mf hys —} TO^ hys + 5ni^\ one gets the second order contributions since the only real diagrams of order 
n are one-loop diagrams in which a counterterm of order n — 1 is inserted. 

We will use the on-shell renormalisation scheme, where the mass and the wave-function counterterms cancel 
all loop contributions to the self-energy such that the pole of the lepton propagator is equal to the physical mass 
of the lepton. Then we obtain the relation 

m (°) = mfy* + y 6m { , n) = p— + ^— . (19) 
' 1 ^ 1 1 + A| 1 + A| y J 

n—l 

For the numerical analysis, let us consider the mSUGRA scenario SPS4; its parameter values are given in 
Table [TJ For this model, the constraint reads (Figure [4]) 

. „ „ . f 0.097, if <&V/ 8 / > 

\Srr6ll\ <\ (20) 
1 RRLLl \0.083, ifM<0 
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Figure 4: Percentage deviation of the electron mass through SUSY loop correction in dependence of real S R 3 R and 
5\\ forSPS4. 



scenario 


x = 0.3 


x = 1 


x = 1.5 


x = 3.0 


for 


1 


Mi = Mi = rriL = rriR 


0.261 
0.234 


0.073 
0.059 


0.050 
0.040 


0.026 
0.023 


^Vii > o 
^Vii < o 


2 


3Mi = Mi = rriL = trr 


0.301 
0.269 


0.083 
0.067 


0.057 
0.045 


0.029 
0.024 


> o 


3 


Mi = M% = 'irriL = rriR 


0.292 
0.235 


0.082 
0.067 


0.057 
0.042 


0.031 
0.027 


> o 

M < o 


4 


Mi = M 2 = 2f - m R 


0.734 
0.702 


0.210 
0.190 


0.142 
0.127 


0.071 
0.064 


> o 

M < o 


5 


3Mi = M 2 = m L = 3m R 


0.731 
0.693 


0.205 
0.179 


0.137 
0.116 


0.067 
0.054 


> o 



Table 2: Different mass scenarios and the corresponding upper bounds for \S R 3 R S]^ L \. m R Lj denotes the average 
right and left-handed slepton mass, respectively, Mi and M 2 the bino and wino masses. In all scenarios, tan j3 = 50 
and sgn(/i) = +1. The upper line is valid for S R 3 R Sj^ L > 0, the lower for Sj^Sj^ < 0. 



As discussed above, the constraint is independent of the overall SUSY mass scale. It does, however, depend on 
the ratio of the various masses, in particular x = ii/mn.. tur^ denotes the average right and left-handed slepton 
mass matrix, respectively. To be sure that this result is no special feature of SPS4, we consider other scenarios in 
Table[2l which differ by ratios of the SUSY breaking masses. The sparticle mass spectrum at the electroweak scale 
for SPS4 corresponds most likely to scenario 2 with x ~ 0. 9 — 1. The bounds are very strong for larger values of 
x and weaken for a small ^-parameter. The upper bounds do not change considerably for values of x larger than 
one and therefore are relatively stable and independent of the parameter space. Electroweak symmetry breaking 
yields a relation between the /x-Parameter, the mass of the Z boson and the two Higgs fields such that in absence 
of any fine tuning u 2 should be within an order of magnitude of M§ (also known as the ii problem). 

We will see in the following section that the dominating RR terms in the flavor violating self-energy r —> e 
cancel in part of the parameter space. That is why so far no upper bound on S RR could be derived |11|J12|. In 
these regions, we can use the constraint on I^r/j^IlI as a restriction on 6 R 3 R . 

3.2. Lepton-flavour violating self-energies 

Lepton flavour violating self-energies can be tan /3-enhanced and, moreover, they also have a non-decoupling 
behaviour. They occur in the renormalisation of the PMNS matrix and lead to a correction of the radiative 
decays lj — > ^7. 
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r v -<SKe £ f L/ -'^e L u T/ -^u e 

e-L TR \ f e L t r \ f e L r R \ f e L 

_ffV__ p J 7 S ffVjXff H^^_p<fW 

X v « X"" X"» 

Figure 5: Dominant diagrams for the neutralino-slepton and chargino-sneutrino loop for the tr — > e^-transition. 

We consider all diagrams with one LFV MI in the slepton propagator and start with tr — ► e^; the dominant 
diagrams are shown in Figure [5] In fact, we can do an approximate diagonalisation of the neutralino mass matrix 
and use the MI approximation in the slepton propagator. Furthermore, we will choose fi to be real. 

The dominant diagrams are proportional to the mass of the tau and sensitive to the LL element (see diagram 
1 to 3 in Fig. (0), while the RR dependence is suppressed, 

s£-« c fx tan/3 m, L m fL 5f L {~M X h (M*, m\ L , m% R ,m% L ) 

+ i^M 2 ) h {M\, fi 2 , m% L , my | . (21) 

This self-energy will contribute to the renormalisation of the PMNS matrix (Section 13. 3p . As in the previous 
section, this contribution potentially leads to an upper bound on 5\f L when a naturalness argument is applied to 
the PMNS element U e 3- Note that in Eq. (|2~Tj) the LFC mass renormalisation is already taken into account. 

The dominant contributions for the opposite helicity transition, tl — > £r, are analogous to the first and second 
diagrams in Fig. [51 

-o cti m phys 

^r L -e R ~ ^ YTA7 Ml M tm ^ ^ 5 * R ^ ' ^ ' ' ) _ H ^ ' ' ' ^ ^ ) ' (22) 

They are sensitive to the RR element; however, the relative minus sign due to the different hypercharges potentially 
leads to cancellations. In this approximation, the RR sensitivity vanishes completely for fi = rrif L and hence no 
upper bounds on S % RR has been derived, as mentioned in the previous section. In Ref. [12] the processes it, — > and 
fx — > e conversion in nuclei are combined. The corresponding one-loop amplitudes suffer from similar cancellations, 
albeit in different regions of the parameter space leading to the constraint S RR < 0. 2. 

Let us now turn to the chargino-sneutrino loops. As the left-handed charged sleptons and sneutrinos form a 
doublet, they have the same SUSY breaking soft mass and therefore the same off-diagonal elements. The neutrino 
is always left-handed so that the chargino loop can only be sensitive to the LL element and a chirality flip of the 
charged lepton is needed. The higgsino component of the chargino couples to the right-handed lepton and the 
wino part to the left-handed lepton. Thus, the self-energy is proportional to the mass of the right-handed lepton, 



3 



= 1^5*5 E E Z^Zrzrm^Bo (m? nJ m|, 

n=l k=l,2 

Ct 2 m 



47T 



M 2 ll tan mv, m Sj 8% L h {rn% ,m% . ,M 2 V) . (23) 



In the second line, we used the MI approximation in the chargino propagator and for the LFV (see last diagram 
in Fig. fl5}). 



3.3. PMNS matrix renormalisation 

Up to now, we have only an upper bound for the matrix element U e 3 and thus for the mixing angle #13; the 
best- fit value is at or close to zero (cf. Eq. (|7J|). It might well be that it vanishes at tree level due to a particular 
symmetry and obtains a non-zero value due to corrections. So we can ask the question if threshold corrections to 
the PMNS matrix could spoil the prediction 13 = 0° at the weak scale. What does it mean for the physics at the 
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high scale if experiment will tell us that #13 does not vanish? As before, we demand the absence of fine-tuning 
for these corrections and therefore require that the SUSY loop contributions do not exceed the value of U e 3, 



<■.', 



< 



U 



phys 



e3 



(24) 



13 
LL- 



Then we can in principle use the smallness of U e 3 to constrain 8 

In an effective field theory approach, the renormalisation of the PMNS matrix is done via rotation matrices 
that diagonalise the mass matrix, which receives contribution from both the tree-level coupling of the fermions to 
Hd and the loop-induced coupling to H u (see Refs. [13,21.28,29] for the quark sector). Lepton flavour violating 
self-energies induce off-diagonal entries in the mass matrix. In order to deal with physical fields, one has to rotate 
them in flavour space to achieve a diagonal mass matrix. 

As a drawback, the effective field theory method is only valid if the masses 
of the supersymmetric particles in the loop are much larger than v — 174 GeV. 
For sleptons and neutralinos this assumption is doubtful, so that we resort to 
the diagrammatic method of Refs. [T6HT8] which does not rely on any hierarchy 
between MgusY and v. In our diagrammatic approach, we consider chargino 
and neutralino loops in the external lepton propagator and resum all tan (3- 
enhanced corrections explicitly. Once again the on-shell scheme is used. The 
loop corrections are finite and the counterterms are defined such that they exactly 
cancel the loop diagrams: 

= u vhys + su^ — u phys 



su 




(25) 



The first-order correction is displayed in the adjoining figure. The counterterm reads 



SUP = Yuf k hys - 



mi 



jk 



'v rrphys 1 
l^ljtj u lk mj 



y> j-rphys 1 
l^tyj U lk mi 



j > I , i.e., heavy particle as external leg; 

j < I , i.e., heavy particle as internal propagator. 



(26) 



As for the mass renormalisation there are no genuine tan /3-enhanced two-loop diagrams. The corrections in 
second order come from one-loop diagrams in which a counterterm of first order is inserted, corresponding to the 
substitution Uf k hys — > Uf^^ + SUfP . In contrast to the resummation of the mass counterterms, these counterterms 
are not directly proportional to the PMNS-element under consideration. The sum of the counterterms has to 
cancel the corrections up to that order, so at the n th order, one gets 
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(27) 



Now we can take the limit n — > 00 and obtain a linear system of equations for the U^ elements (A; = 1, 2, 3): 
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(28a) 
(28b) 
(28c) 



In the MSSM, we have E 



E*° 



E^ . As shown above, E Ti? _ ei is sensitive to 5j^ L and so is SU e 3. We aim 
to avoid accidental cancellations and set all off-diagonal elements to zero except for Jj^JE In this case we can 
explicitly solve the linear system of equations 
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(29) 



2 In 1261 we study the case with nonvanishing S]^ R - 
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Figure 6: \8U 13 \ /U^ ys in percent as a function of 5j^ L and tan/3 for 8 13 = 3° (left) and 9i 3 — 1° (right). 



By means of Eq. (|24p . we can in principle derive upper bounds for 5]j L . As shown in Figs. [6] they strongly depend 
on tan/3 and the assumed value for U^ ys . 

After three years of running, the DOUBLE CHOOZ experiment will be sensitive to 6*13 = 3°, which corresponds 
to U e 3 = 0.05. A future neutrino factory may probe #13 down to #13 = 0.6° [3D]. In general, even with future 
experimental facilities, we can conclude that the corrections from SUSY loops to the small element U e 3 stay 
unobservably small. This means at the same time that if some experiment measures #13 ^ 0, this will not be 
compatible with tri-bimaximal mixing at the high scale and moderate sparticle masses, since SUSY threshold 
corrections cannot account for such an effect: Even for large tan/3 the already existing constraints on 8)^ L from 
t — y s'y are stronger assuming reasonable SUSY masses. However, since r — > decouples, our method leads to 
a sharper bound for very large SUSY masses, especially with 6*13 = 1° and large tan/3. 

3.4. Counterterms in the flavour basis and charged Higgs couplings 

Neutrinos are both produced and detected as flavour eigenstates. In order to have flavour diagonal W couplings, 
however, it is necessary to introduce counterterms, SVij, which cancel the LFV loops. By doing this you perform 
a renormalisation of the unit matrix. In an effective field theory approach this is achieved via a wave function 
renormalisation by rotating the lepton fields leading to a diagonal mass matrix and physical fields. This rotation 
of the fields induce LFV in the charged Higgs coupling to lepton and neutrino and the same is true for the 
counterterms in the diagrammatic approach. 

The first-order correction is displayed in the figure above; the flavour-diagonal vertices do not get any coun- 
terterms, since the external loops are already included in the mass renormalisation. We obtain 
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(30) 



(0) 



These counterterms 



You can translate this to the mass eigenstate basis used in Eq. P5|) via SUik = ^i*jUj k 
induce LFV in the charged Higgs coupling to lepton and neutrino, due to the different helicity structure of the 
Higgs and W coupling and the different lepton masses. The H + ev T vertex can be of particular importance, since 
it is possible to pick up terms with a tau Yukawa coupling. As discussed before, this coupling is enhanced in the 
large tan/3 regime and can partly compensate the loop suppression factor. 

The chargino contributions from the counterterm and the LFV loop cancel in the charged Higgs coupling as 
the chargino loop is exactly proportional to the mass of the right handed lepton. Therefore only the neutralino 
contributions remain. 
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The charged Higgs coupling to electrons reads 

V2M W ^ V m? hys / \/2A% M V m? hys 1 + A e 1 + Ar J V ' 



^ert—TL j 


«52 


m? hys J 




Z-'eR—fiL ] 


1 _ ^2 







tan ^ H 0) ^ + = ^ tan H"3S^ + fee 1 • (31b) 

^ " 1 + A e ^ V 1 + A r Alr ) ' (31C) 

where Eg V = tan/3 A^ R . We see that the counterterms are suppressed with the electron mass. As the 

lepton mass cancels out in H eR _g iL /m^, the LFV loop contributions with v T and differ by a factor Tot/to^. 
Similarly, we obtain for the coupling to muons, 

tT *£ = tan/3 ( m (o) §v + m (0)h^) = _^2_ p f_ m^»^ ^^-) , (32a) 

g+ _ *92 < ys / tan/3 \ 

^ " V2A% 1 + A^ ^ m P^ 1 + A T LR J ^ 2h > 

< - 7& " <S?) = tan ^ " ^Efe:) ■ (32c) 

While the first term is similar to the couplings to electrons, the counterterm dominates over the loop contribution 
if there is an electron neutrino in the final state. 



Finally, for the r coupling one finds 

z phys 

y/2M w 1 + A T 



ir£ = -J^- f^-j- tan /3 , (33a) 
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tan /3 


1 TR — /i£, 


m? hys J 


V2M W 
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771? hyS J 


y/2M w 


<1 + A T 



,rf; = tan/3 f m (?)*y 3 i - ) = ta / £^ _ mg^ S \ 

™ c V2M W ^ V m? hys / V2A/ W \ 1 + A T 1 + A e / 1 y 

The results of Eqs. ([5T ]) -(|gg ]l are given in Eqs. (92-95) of Ref. [S3] for the decoupling limit M SU SY > v. In 
Appendix C of Ref. [23J an iterative procedure (analogous to the one in Ref. [H]) has been outlined which achieves 
the all-order resummation of the tan /3-enhanced higher-order corrections. Eqs. ([3T|) - (|33|) comprise analytical 
formulae for the limits to which this iterative procedure converges. 

The tan /3-enhanced lepton flavour violating Higgs couplings can become important in the leptonic decay of 
charged Kaons, K — > li>, where they potentially induce lepton non-universality. Then the current experimental 
data and our fine-tuning argument together constrain the various terms in Eqs. (|31j) . as they contribute to the 
electron self-energy as well. In particular, if the second term in Eq. (|31c[) had a significant effect in the ratio 
Rk = V(K — > ev)/T{K — > (iv), as was assumed in Ref. [24], & e LR would give a large contribution to the electron 
mass [3T]. (The value A^ = A e LR = 1(T 4 [24] corresponds to 

^rr&lJl ~ 2 in the SPS4 scenario and thus gives 
a more than 2000% correction to the electron mass.) While in the improved analysis [32] the contribution of 
^to the electron mass was not considered, their scanned values of 5j^ L RR are in agreement with the fine- 
tuning argument. The scan respects ^ix,^^ < 0. 25, in marginal agreement with our results of Sects. 3.1 and 
3.5. The NA62 experiment at CERN aims to reduce the error of Rk from 1. 3% to 0. 3%. This prospective error 
is used in Ref. [33] to derive large, phenomenologically interesting values for 5j? L and 

3.5. Anomalous Magnetic Moment of the Electron 

The anomalous magnetic moment of the electron plays a central role in quantum electrodynamics. The precise 
measurements provide the best source of the fine structure constant a em if one assumes the validity of QED [31] . 
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Conversely, one can use a value of a em from a (less precise) measurement and insert it into the theory prediction 
for a e to probe new physics in the latter quantity. The most recent calculation yields |35| 



1 159 652 182.79(7.71) x 10 



-12 



(34) 



where the largest uncertainty comes from the second-best measurement of a cm which is a"^ = 137. 03599884(91) 
from a Rubidium atom experiment [36J. 

Supersymmetric contributions to the magnetic moment are usually small, due to the smallness of the electron 
Yukawa coupling and the SUSY mass suppression. However, multiple flavour changes, resulting in a LFC loop, 
insert the r Yukawa coupling, which strongly enhances the amplitude. As a result, supersymmetric contributions 
can be as large as 0(1O~ 12 ), comparable to the weak or hadronic contributions [35]. The amplitude can exceed 
a 3cr deviation of the theoretical mean value, which enables us to constrain the LFV parameters 5}^ L and S]^ R . 
In Ref. [TT] the magnetic and electric dipole moments and di of the charged lepton li were calculated in the 
MSSM, considering flavour-conserving and flavour-violating contributions within the mass insertion approxima- 
tion. The authors found that the naive mass scaling can be overcome with double mass insertions. However, in 
their phenomenological analysis to constrain the flavour- violating parameters <5^ Y , they only used a M and the 
experimental bounds on and d e but did not consider a e . Our consideration of a e adds a novel aspect to 
the phenomenological study of LFV parameters in the MSSM and complements the analysis of Ref. [TT] in this 
respect. 

The supersymmetric contributions to the anomalous magnetic moment a e are generated by chargino and neu- 
tralino loops, where the photon couples to any charged particle in the loop. The full analytic result can be found 
in Ref. [37]. Here, we will neglect the terms which are both proportional to the electron mass and not (potentially) 
tan /3-enhanced and are therefore left with 
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The loop functions are listed in Eq. (|80|) and the couplings read [38] 
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(36c) 
(36d) 



The mixing matrices are defined in Appendix [Bj Note that they are 6x6 matrices, in order to allow for flavour 
changes in the loop. 

The dependence on tan/3 in Eqs. (|35[) is hidden in the mixing matrices. In principle, tan/3 comes from a 
chirality flip on the selectron line and in the chargino case from the combination of vacuum expectation value v u 
and the Yukawa coupling, y e v u = m e tan/3. We can, however, simplify the expressions significantly as follows: 
We assume a universal SUSY mass, real parameters and the same signs for Mi and M 2 [39J, then expand a e in 
powers of M^/MgusY or l/tan/3. Then we obtain 
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(37) 



The result is again finite. The 1/M| USY dependence reflects the decoupling behaviour of supersymmetry. Fur- 
thermore, we note that a large value for tan/3 can dilute the l/-A-iJ USY suppression. The numerical results are 
computed with the exact formula in Eqs. (|35[) . 
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Figure 7: Supersymmetric contributions to a e as a function of 5jj L for S R 3 R for (a) scenario 5 (from steep to 
level: S R 3 R = 0.6 (green); 0.4 (blue); 0.2 (red)); (b) scenario 2 (S R 3 R = 0.6 (green); 0.2 (red)) of Tab. with 
Msusy = 500, tan j3 — 50, and /i = Msusy- The light, medium, and dark grey regions correspond to the 
theoretical la, 2a, and 3<r regions, respectively. In (b), the dashed curve shows the result without the mass 
correction. 



So far, the Yukawa couplings are unrenormalised; the inclusion of the mass renormalisation amounts to a loop 
contribution to a e which approximately grows as tan 2 (3 |17) . Diagonalising the mixing matrices perturbatively, 
one finds a linear dependence on the Yukawa coupling of the remaining second terms of Eqs. (|35p . In this way 
we find an easy expression, which takes the corrections into account by a global factor, 



„SUSY, \L , SUSY,A e 
(*„ ~~r (-1, 



.SUSY, 1L 



1 



1 + A f 



(38) 



where af USY ' 1L = a}° + a* , as discussed in Ref. [T7] , 

For the numerical analysis, we only allow 5}j L and S R 3 R to be non-zero such that they are the only source of 



flavour violation. The theoretical uncertainty in Eq. (|34|l is taken as la deviation and we require that the SUSY 
contribution to a e is less than 3a. 

We show the results for our scenarios 2 and 5 (see Table [5]) in Fig. [7] As 8 R 3 R increases, the bound on 5j^ L 
becomes stronger and vice versa. The bound strongly depends on the SUSY mass. Since a e decouples for large 
SUSY masses, the bounds become very loose for Msusy ^ 500 GeV. On the other hand, small SUSY masses lead 
to complex slepton masses, resulting in a lower bound on the SUSY mass. For this reason, the upper bounds 
on a nd S R 3 R are limited by the SUSY mass constraints. We find |<5]^ ■ < 0. 1 for AfsusY £ 500 GeV, 
coinciding with our non-decoupling bound in Eq. (|20[) . 



3.6. The radiative decay lj — >• ^7 

Since their SM branching ratios are tiny, supersymmetric contributions to lepton flavour violating decays ij — > lj"/ 
can be sizable and vastly dominate over the SM values. As indicated above, these decays currently give the 
best constraints on the left-left (LL) and left-right (LR) lepton flavour violating parameters. At one-loop level 
and within MIA, ij — » Ijj has for example extensively been studied in Ref. [11], constraining e.g. the mSUGRA 
parameters Mi and tcir. In this section, we compute the supersymmetric contributions to ij — > Ijj, including 
both the mass renormalisation and the two-loop contributions coming from flavour-violating loops. The current 
upper bounds for the branching ratios are listed in Table [3l 
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Let us briefly summarise the formalism. Three SUSY diagrams contribute to the amplitude of lj — > l^, 
corresponding to the coupling of the photon to lj, li, and the charged particle in the loop. The off shell amplitude 
can be written as [41] 

iM = iee»*u t (p - q) [q 2 iMi P L + AfP R ) + m h %o ^ \A\P h + A R P R )] Uj (p) , (39) 

where e* is the photon polarisation vector. If the photon is on shell, the first part of the off-shell amplitude 
vanishes. 

The coefficients A contain chargino and neutralino contributions, 

A L,R = A (x°)L,R + A (x±)L, Rj , = 1j2j (40) 



so A L is given by the sum of [32] 
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(41) 
(42) 



with the couplings given in Eqs. (f3"6")l . We get A R by simply interchanging L -H* R. 
Finally, the decay rate is given by 

T(l s -+ la) = ^% (\A L 2 \ 2 + \A*f) . (43) 

Both the flavor-conserving transition ij — > Ij j and the flavour-changing self-energies are tan /3-enhanced. For this 
reason, we do not only consider the effect of the mass renormalisation but also include the two-loop contributions. 
Because of the double tan /3 enhancement they can compete with the first non- vanishing contribution. As for the 
corresponding counterterms, mass counterterms have to be inserted. In addition, wave-function renormalisation 
counterterms play a role as the above-quoted result for lj — > 1^ presumes an expansion in the external momenta of 
the lepton. Therefore, to be consistent, the counterterm has to be given in higher order of the external momentum. 
However, only the mass counterterm will be tan /3-enhanced because of the chirality flip involved. Corresponding 
diagrams are shown in Figs. [5] 

The wave-function and mass counterterms are given by: 



1°l = ( 1 + *Sl L l° R = ( 1 + ^Sl R )l R , to° = mi + Smi , (44) 



\SI L ) II , l° R = (l + \c 

where the fields and masses with a superscript are the unrenormalised fields. In order to identify the countert- 
erms, one first considers the kinetic and the mass term of the Lagrangian. The one-loop self-energy of the lepton 
can be divided into a scalar and a vector-type part, where the latter can further be divided in a left-left and a 
right-right transition, 

*E { (p) - iX? L _ lR (p)+ipZ lL - lL (p)P L +ipZ lR - lR (p)PR . (45) 

Now we demand that the additional terms in the mass Lagrangian cancel the scalar-part of the one-loop self- 
energy whereas the additional terms in the wave- function Lagrangian cancel the vector- type part. Therefore the 
counterterms have to fulfil the following conditions: 

Sl L = -Y,i L _i L (p 2 = mf), 5l R = -£i R _ iR (p 2 = mf) , (46) 

5 mi = Z lL - ln (p 2 = mf) - ?±(6l L + Sl R ) . (47) 
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Table 3: Current upper bounds for BR(lj — > hj), j > i \40j . 
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Figure 8: Two-loop contributions to r — > e7 from (a) chargino and neutralino loops; (b) from the counterterm of 
the t propagator. 
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Table 4: Upper bounds on \S]^ L \, Iflj^l and \6j? L \ for the mSUGRA scenarios of Table [TJ from BRll^ — > ^7) 
including mass renormalisation and two loop contributions. 



To give an explicit expression for the counterterms, we expand the self-energies up to the quadratic order in the 
external momentum and then compute the two parts of the one-loop self-energy. In the series, the even and 
odd orders contribute to the scalar and the vector type part, respectively. The chargino contribution to the 
counterterm is then given by 

Sm = i^2E r, ' 5 '» r ^» f--Bo + mfm & C (m|,m|,m|)- -mfm Xfc £> 2 (m|,m|,m|,m|) 

k ^ 

= 1^2 E r ^.** r ^4 aj(ro *' m * >m|) 

k 

k 

where d = 4 — 2e. The wave-function counterterms also induce an additional lepton photon vertex, 81 l IlI^IlA^. 

Now we can compute the various diagrams (cf. Fig. [5]). Up to second order in the momentum p all contributions 
indeed cancel each other. For the chargino two-loop contribution to t — > ej, we obtain 
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where 



^ax = - (\Ntuc\ 2 + |iV 3 V| 2 ) Jz F i N M ~^ Rc [NsAxNflx] F 2 N (x AX ) , 

/ 7 ?\ 1 2m~± 

%x = (|C 3 V | + |C 3 V I ) j^Ff ( XAX ) + -3^ Re [CJ^C©*] i^ 2 c (^x) 

and the couplings N^J^ and Cjj^jr are defined in Eqs. (|36p . 

For the numerical analysis, we first consider the mSUGRA scenarios listed in Table [1] as well as the scenarios of 
Table [21 The ti parameter at M ew is determined with Isajet |43H45| . Note that the different bounds for 8 l [ L in the 
literature can differ due to their dependence on the chosen point in the SUSY parameter space, see, e.g., Refs. [411 
|41^|41)1|4"T] . Table U] summarises the bounds on |<5]^|, \^ll\ an< ^ |^£l| f° r these scenarios; they include both tan/3- 
enhanced corrections to lj ^7, namely the mass renormalisation and two loops contributions. Interestingly, the 
two corrections tend to cancel each other: As illustrated for SPS4 in Table \5[ the mass renormalisation tightens 
the bound, whereas the two loops effects increases them again. Thus, the effect is generally smaller than 1 %, 
particularly for the small tan/3 scenarios. For large tan/3 (as is the case in SPS4), however, the deviation can 
reach 6 %. Without the inclusion of our two corrections we recover the results already found in Refs. |11[[P2] after 
taking into account that the experimental upper bounds have changed a bit. 

Let us therefore study the scenarios of Table [2] with Msusy = 300 GeV and tan/3 = 50. The branching ratios 
of lj — > kj in the scenarios 1, 2 and 5 of table [2] are shown in Figure [9l We see that scenario 5 gives the strongest 
constraint on 8 13 L . In addition, the corrections discussed here have the biggest effect in this scenario. The upper 
bound on Sj^ again depends on the SUSY mass. The branching ratio decouples for large SUSY masses so that 
the upper bounds weakens for increasing Msusy (Fig. fTO]) . 

As already noted in Sec. 13.31 the corrections from supersymmetric loops cannot reasonably push U e 3 into the 
reach of the DOUBLE CHOOZ experiment without violating the bound from r — > ej: E.g. for sparticle masses of 
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Figure 9: The branching ratio as a function of 6}^ L , 8 2 ^ L and 6j? L with corrections in the scenarios 1, 2, 5 (bottom 
to top) at Msusy = 300 GeV and tan f3 = 50. 
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Table 5: Upper bounds on \Sj^ L \, \S}^ L \ and \Sj^ L \ for SPS4 without any corrections, with mass renormalisation and 
taking into account both mass renormalisation and two loops contribution. In parenthesis: deviation compared 
to the tree level bound in percent. 
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500 GeV we find \SU e 3 \ < 10 -3 corresponding to a correction to the mixing angle #13 of at most 0.06 degrees. That 
is, if the DOUBLE CHOOZ experiment measures U e s 7^ 0, one will not be able to ascribe this result to the SUSY 
breaking sector. Stated positively, U e 3 > 10~ 3 will imply that at low energies the flavour symmetries imposed on 
the Yukawa sector to motivate tri-bimaximal mixing are violated. If the same consideration is made for the most 
optimistic reach #13 < 0. 6° of a future neutrino factory (the quoted bound corresponds to the best value of the 
CP phase in the PMNS matrix), the threshold corrections become relevant only for sparticle masses well above 
1500 GeV. While the considered effects in both r — > and U e s involve the product Sj^ L 5}^ L , the qualitative result 
is equally valid, if the needed flavour and chirality violations are triggered by 6^ L or other combinations of the 

While we have only considered loops with a single flavour-changing d^y m our discussion of lj — > ^7 decays, 
contributions proportional to Sj^ R Sj} L can be relevant for [i — > ej. For recent analyses including this effect we 
refer to Refs. |48U49| . If both tan/3 is large and the charged-Higgs-boson mass is small, further two- loop effects 
involving a virtual Higgs boson can be relevant [SHIHH ■ These effects are suppressed by one power of tan [3 with 
respect to the two-loop corrections included by us, but do not vanish for A/susy 3> Mh+,v. In Ref. [23] effective 
lepton-slepton-gaugino vertices reproducing the chirally enhanced corrections in the leading order of w 2 /A/| USY 
have been derived and applied to electric dipole moments, cf. the overview on previous work at the beginning of 
Seel 



4. Renormalisation group equation with SUSY seesaw mechanism 

In the previous section, we derived bounds on the off-diagonal elements of the slepton mass matrix, parametrised 
by S^y These are a priori free parameters in the MSSM; they are set once we know how supersymmetry is 
broken. We saw, however, that these elements are well-constrained and this result generally applies to the soft 
terms. Therefore one usually assumes universality of the supersymmetry breaking terms at a high scale, e.g. 
Mgut = 2 • 10 16 GeV where the SM gauge couplings converge. Then the renormalisation group equation (RGE) 
running induces non-vanishing S^ Y a t the electroweak scale. Clearly, the size of 8 l ^ Y IS model-dependent. 
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In this section, we will consider two GUT scenarios based on the gauge group SO(10), which generically includes 
right-handed neutrinos. The breaking of SO(10) around Mqtjt generates heavy Majorana masses for these right- 
handed neutrinos. After the electroweak symmetry breaking, the left-handed neutrinos receive small Majorana 
masses via the seesaw mechanism. 

4.1. Neutrino Yukawa Couplings and Grand Unification 

The seesaw mechanism naturally explains tiny neutrino masses. As already discussed in the Introduction, the 
right-handed neutrinos are singlets under the standard model group. Then we can write down an explicit mass 
term, (Mii)ijvfVj (see Eqs. ([1])). Now if the entries (Mr)^ are much larger than the electroweak scale, we can 
integrate out the heavy neutrino fields at their mass scales. Below the scale of the lightest state the Yukawa 
couplings are then given by 

Weff = Wmssm + \ {Y V LH U ) T M- 1 (Y V LH U ) . (50) 

After electroweak symmetry breaking, W c s leads to the following effective mass matrix for the light neutrinos: 

M v = -Yj M^ x Y v vI = -kvI . (51) 

Since the light neutrinos cannot be heavier than 1 eV and the mass scale of the atmospheric oscillations is of 
order 0.1 eV, the Majorana mass scale is around 10 14 GeV. 

In the MSSM, it is convenient to choose both the Yukawa coupling matrix of the charged leptons and the 
Majorana mass matrix of the right-handed neutrinos diagonal. In this basis, Ai v is diagonalised by the PMNS 
matrix, 

^pmns-^^pmns = diag (m^,^,^ ) = -T> K vl . (52) 

By means of Eqs. ([ST]) and (|52|). Y v can be expressed as [52] 

Y„ = V^jj R Ul MNS , V M = diag (M Rl , Mr 2 , M Rs ) , (53) 

with an arbitrary orthogonal matrix R. Thus Y v depends both on the measurable parameters, contained in the 
diagonal matrix T> K and E/pmnS; an d the model-dependent parameters, namely three Majorana masses and three 
mixing parameters. In the MSSM, these are completely free parameters. 

The seesaw mechanism is automatic in grand-unified models with broken U(l) fl _r symmetry [53]. (B and L 
denote baryon and lepton number, respectively.) In SO(10), the SM fermions of each generation are unified in one 
matter representation, together with the singlet neutrino (541155] . No further fermionic multiplets are needed. An 
additional Higgs field acquires a vev, breaking the SO(10) subgroup SU(2)^ x U(1) B _ £ to hypercharge, U(l) y . 
At the same time, Majorana masses for the SM singlets are generated. As indicated in Eqs. ([1]), the up-quarks 
and neutrinos couple to the same Higgs fields H u so that the Yukawa matrices Y v and Y u are related. The actual 
form of this relation is model-dependent; however, there are two extreme cases |46[I56|. 

1. Minimal (CKM) case: The mixing in Y v is small and 

y v = y u = y c T KM v u u CKM (54) 

holds at the GUT scale. This case refers to minimal SO(10) scenarios with small mixing angles for the 
Dirac mass matrices. The large leptonic mixing angles are a consequence of the interplay of Y v and Mr in 
the seesaw mechanism. 

For normal-hierarchical neutrino masses, i.e. m Ul -C m U2 ~ yj Am^i *C m Vz ~ \J Am^ and the MNS matrix 
being close to its tri-bimaximal form, the masses of the right handed neutrinos are given by 

1 

Mr.i ~ Tn 77? — j Mr 2 sa 2 h 3 , 

3m,2 2777,2 



6m, 



(55) 
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Figure 11: Illustration of the procedure used to solve the RGE 



2. Maximal (PMNS) case: Large mixing in Y v is achieved in models with 



t 

PMNS ' 



(56) 



This scenario is the analogon to the quark case: the mixing matrix arises in the Dirac couplings, with 
the Majorana matrix being diagonal. Note that Y v is not symmetric any more and this relation is indeed 
realised in models with lopsided mass matrices. In this case, the masses for the right handed neutrinos are 
simply 



M Rl = M R2 = ^, M Rs = —L. 

m,M to„ to,,, 



(57) 



In terms of the parametrisation (|53j) , the second case corresponds to R = 1. Then the mixing in Y v is determined 
by the PMNS matrix. By contrast, small CKM-like mixing in Y v (case 1) requires a non-trivial structure of R. 

Clearly, these two cases are special; however, they provide two well-motivated but distinct scenarios. A more 
detailed introduction to these two cases are given in |46[f56] . Note that the authors use the LR convention for the 
neutrino Yukawa coupling so that their equations differ from ours by the substitution Y v <-> Yj . 



4.2. Renormalisation-group Analysis 

We list the renormalisation group equations (RGE) of the MSSM [57H59| and the MSSM with right-handed 
neutrinos [52.60,61] in Appendix [Cl The right-handed neutrinos are singlets under the SM gauge group so that 
they do not change the RGE for gauge couplings and gaugino masses. 

The procedure of solving the RGE is schematically depicted in Figure [TTJ With the experimental values of the 
indicated parameters at the scale Mz, we evaluate the gauge and Yukawa couplings at the various mass scales. 
The three heavy neutrinos are included step by step. At the GUT scale, Mgut = 2 • 10 16 GeV, we assume 
universality of the supersymmetry breaking soft parameters, 
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Solving the RGE in leading order, one gets for the LFV off-diagonal elements [52,62 

3mg + 41 / vt , ( Mgut \ 



A-i)--^^^m^y^, Mil (59) 

Am:),, -0. i^j, (60) 
(^^-A^felnf^V^ - Mi- (61) 



/« - 87r 2 « y Mil . y 

(Recall that m| and m| contribute to (m 2 )z,z, and (m 2 L )nR as shown in Eqs. (j87b[) and (|87c[l .) The size of LFV 
depends essentially on the structure and magnitude of the neutrino Yukawa coupling, the scale of the right handed 
neutrinos as well as on the SUSY breaking parameters mo and Aq. The only source of flavour violation stems from 
Y v . According to Eq. (|60p . off-diagonal RR elements 5 RR are not generated at leading order. The off-diagonal 
elements (.A;)^ , ■ are related to S l ^ R , as can be read off from the slepton mass matrix listed in Appendix [Bj They 
are proportional to Aq and suppressed by a Yukawa coupling. So in general, the generated 8 % l R at the weak scale 
are negligibly small compared to the generated 8^ L elements (see, however, Section EO| . 

In view of the slepton mass matrix, the RGE are usually solved by integrating out the right-handed neutrinos at 
one scale Mr ~ O(10 13 — 10 14 ) GeV. In most GUT models, however, the heavy neutrinos are strongly hierarchical 
so that these degrees of freedom should be integrated out successively. As a result, we have a number of effective 
field theories below the GUT scale; the details are listed in Appendix [Cl The running of the mixing angle can 
change significantly with three non-degenerate heavy neutrinos and cannot be reproduced if all heavy neutrinos 
are integrated out at a common scale M- mt [6T1I63) . 

Our input values are the gauge couplings, the masses of the leptons and quarks at the electroweak scale, the 
neutrino mass differences Am 2 tm and Am 2 ol as well as the PMNS matrix. We will assume the normal hierarchy 
for the masses of the light neutrinos, 



m\ + Am s 2 ol , m„ 3 =^J m l i + Am 2 tm , (62) 

so we are left with the mass of the lightest neutrino m Vl , the three masses of the Majorana neutrinos Mr 4 , the 
factor tan j3 and the uncertainties in the mixing angles of the PMNS matrix as free parameters. For numerical 
results, we will choose m Vl rj 0(1 -3 ) eV. Note that the heavy neutrino masses Mr are already fixed by Eqs. ([55]) 
and (|57|) . 

/(4.0-I0 9 GeV, 4.0-I0 9 GeV, 5. 9 • I0 14 GeV) , PMNS case; 

Mo, , Mb,, Mb, I — < ; „ ' (63) 

V 1 2 3 ; \ (2. • I0 6 GeV, 3. 9 • I0 11 GeV, 7. 4 • f 15 GeV) , CKM case. V ; 

In addition, the soft SUSY breaking terms Aq, mo, and mi/2 as well as tan/3 are free parameters. 

4.3. Numerical Results 

We start by considering the Am? entries in Eq. ([56| for the two scenarios. In the PMNS case, we get from 
Eq. K55J) 

/ a 2n _ 3mg + Al ( 2 M GUT 2 Mqut \ . — 2^ 

( Am i)i2 ~ ^2— (Vt U <*U& In + y c U e2 U^ In j + 0(y u ), 



3ml + 


.12 
^0 


8tt 2 




3ml + 


42 
^0 


8tt 2 




3ml + 


J2 
^0 


8tt 2 



(A-!),, - -^i^^C/.sln^HI + 0{y l). (64) 



Mr. 



In these equations, we replaced by U since we already know that SUSY corrections do not spoil a possible 
symmetry at the high scale (Sec. 13. 3|) . Hence, we can neglect the small difference between JJ phys and U^°\ In the 
following, we will distinguish between two different input values at Mgtjt, $13 = 0° and #13 = 3°. 



3 As we assume universality of the SUSY breaking parameters at Mqtjt, we do not receive additional terms due to the coupling to 
coloured Higgs fields as in, e.g., Refs. |33II60I . 
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Table 6: Results for the generated off-diagonal elements 5 % l L at Mz in the PMNS case for the different mSUGRA 
scenarios. We assume 612 = 33° and 6*23 = 45° according to the tri-bimaximal scenario and consider the two cases 
#13 = 0° (top) and 6 13 = 3° at M z (bottom). 



In the 12 and 13 elements, the large top Yukawa coupling compensates the suppression by U e s, i.e., #13. For 
013 = 3°, which is the sensitivity of the DOUBLE CHOOZ experiment, the top contribution dominates, 

ylUezU^hx^ y 2 U e3 U T3 ln^ 

— tt^- w 13000 • U e3 « 650 , — jrr^- pa 12000 • U e3 ps 600 ; (65) 

y c U e 2Un2ln Mft2 y c U e2U T 2n\ Mr ^ 

however, for much smaller angles, the contribution of the second generation needs to be taken into account 
and leading-order RGE are not a good approximation anymore. Note that the dominant contribution to 5j~ L is 
independent of the unknown mixing angle #13. 

We can perform the analogous analysis for the CKM case. Then the generated 8 l [ L are one or two orders of 
magnitudes smaller for this case, simply because the small CKM mixing angles replace the large PMNS ones in 
Y v (see Eq. (HID). 

In the following, we will analyse how large the LFV off-diagonal elements 6^ L can get at the electroweak scale 
due to renormalisation group running and study their sensitivity on 9\ 3 . 

PMNS case. The main contribution stems from the running between Mqut and Mr 3 , where the dominant 
entry of Y v is of the same order as the top Yukawa coupling. Below Mr 3 the entries of the remaining neutrino 
Yukawa coupling is much smaller such that the result is only weakly dependent on the two lighter Majorana 
masses. 

Table [6] lists the results for two different neutrino mixing angles; we generally obtain \Sj^ L \ < |<5ii| < riz-r ^ s 
expected from the leading-order RGE the sizes of the (1,2) and (1,3) elements increase by two or three orders of 
magnitude for a sizable U e3 (Mz) = 0. 05 element, compared to the case with U e3 = 0. For 6*13 = the estimate 
from the leading-order solution does no longer coincide with the exact numerical solution. While for 6*13 = 3° the 
relation 



(Am|) (Am|) 



13 



1 holds, the latter estimate is no longer valid for 9\ 3 = 0. 



Now we use b x L \ in order to derive an upper bound on 9\ 3 for the different mSUGRA scenarios and obtain 

1 013 1 < (0. 25°, 0. 42°, 1. 1°, 2. 2°, 0. 30°, 0. 5°, 1. 2°) (66) 

for the respective scenarios. The element 6j? L is far less sensitive to 6*13: even for of 6*13 = 3°, it is at least one 
order of magnitude below the current experimental bounds. As discussed above, Sj^ L is not sensitive to 0i3. In 
the SPSla and SPS4 scenarios, however, it is above the experimental bound, whereas it is well below the limit in 
SPS2, SPS3 and B. Hence, some region of the parameter space can be excluded by the element 5 2 L \. 



CKM case. The neutrino Yukawa matrix contains an C(l)-entry only above the scale Mr 3 . Thus 11011- vanishing 
8 % l L are basically generated in the interval [Mr 3 ,Mgut]- 

The results are shown in Table [71 As expected, the values for the various 8 % l L are small, due to the small CKM 
mixing angles. They are well below the experimental bound so that we cannot exclude parts of the parameter 
space in this case at all. 

In summary, we have seen that LFV processes offer a window to look at the structure of SUSY GUT scenarios. 
The results of the various cases considered in this paper (see Tab. 01 [S] and [7J are compared in Figure The 
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Table 7: RGE induced off-diagonal elements at Mz 8 % l L in the CKM case for the different mSUGRA scenarios. 
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Figure 12: Summary of the results, as shown in Tables^® and\?\ The red circle denotes the experimental bound, 
the blue square the PMNS case with 813 — 0°, the green diamond the PMNS case with #13 = 3°, and the orange 
triangle the CKM case. 



decay r — > ^7 can exclude the PMNS case through 5j? L for some SUSY mass spectra, irrespective of U e s. In 
addition, #13 is bounded by 5}? L ; the PMNS case allows for small values of #13 only. A more precise measurement 
of U e 3 will make it possible to disfavour or even to exclude models. In the CKM case, /j, — > and r — > ^7 should 
be observable in the near future [M]. Here, the tan /3-enhanced corrections should also be included. 

Most of these conclusions hold irrespective of the GUT correction terms which will be discussed in the following 
subsection. These do not affect the LL sector; however the LR sector will be modified significantly. Furthermore, 
note that any observation of r — > e-f calls for additional sources of LFV; as discussed above, the needed 5j^ L 
cannot be generated if we start with universal boundary conditions at Mqut- 



4.4. Effects from Fermion Mass Corrections 

In grand-unified models, the Yukawa couplings arise from few basic couplings, relating the couplings of the SM 
fields. In particular, minimal GUT models predict the unification of the down quark and charged lepton masses. 
While those of the bottom quark and the tau are in remarkable agreement at Mqut, the relation is violated for 
the first and second generation. The failure for the lighter generations, however, is naturally explained by the 
presence of higher dimensional operators due to physics at the Planck scale that induces corrections of the order 
Mgut/Mp! [55]. 

These nonrenormalisable operators do not only help to get realistic fermion mass relations, they also cure 
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another problem of SUSY SU(5) models: the too large proton decay rate stemming from couplings of the color 
triplet Higgs field (>(> 7X3 • The consequences of the higher-dimensional operators on flavour physics observables, 
however, were for a long time neglected. First studies were carried out in Refs. |71[I72| with vanishing neutrino 
masses. In Ref. [72] the RGE of the RR slepton sector involves the CKM matrix together with an additional 
rotation matrix. Ref. [7T] discusses an SO(10) SUSY GUT model with nonminimal Yukawa interaction. Massive 
neutrinos via the seesaw mechanism were first analysed within SUSY SU(5) in Ref. [73] ■ The authors include 
one higher-dimensional operator whose main effect is described by two mixing angles that parametrise rotations 
between the down-type quarks and charged leptons in the first and second generations. In Ref. [74| the correlation 
between the flavour-violating mass insertions of the squark and slepton sectors in a SUSY SU(5) scenario are 
studied, including the corrections to fix the quark-lepton mass relations. The authors of Ref. [75] study an 
SO(10) SUSY GUT model in which the large atmospheric mixing angle can induce — sr transitions. They 
parametrise the effect of higher-dimensional operators on flavour transitions in terms of a mixing angle and a 
(CP- violating) phase and place tight constraints on these parameters from a simultaneous study of K—K mixing, 
Bd—Bd mixing and B s —B s mixing. A very detailed theoretical analysis generalising the approach of Ref. [73] has 
recently been performed in Ref. [75] , for a compact summary see Ref. [77] . These papers contain a complete list of 
RGEs for SUSY SU(5) including nonrenormalisable operators for all three types of the seesaw mechanism. This 
setup drastically increases the number of free parameters which show up in several diagonalisation matrices. Even 
with flavour-blind and field-type-independent mediation of SUSY breaking, the higher-dimensional operators give 
rise to tree-level flavour-violating entries in the sfermion mass matrices. Their effective trilinear couplings are no 
longer aligned with the effective Yukawa couplings. Since the A-terms contribute to the slepton mass matrices 
already at tree level, these misalignments are potentially very dangerous. The authors of Ref. |76| study special 
types of Kahler potentials and superpotentials in which such terms can be avoided. With some approximations 
they recover the parametrisation with mixing angles between the first and second generation adopted in Ref. |73| . 
We will adopt a similar approach explained below. A comprehensive phenomenological analysis with the RGE of 
Ref. [7B] and its very general diagonalisation matrices has not been done yet. We will make a simplified ansatz: 
Instead of using the most general setup, we concentrate on the lepton sector and parametrise the effect of higher 
dimensional operators as a rotation between the first and second generation without the inclusion of any phases 
(which are not probed by current experiments). Further we only use the RGE of the MSSM and focus on the 
effect in the trilinear terms. 

If we denote the renormalisable and the higher-dimensional couplings as Ygjjt and Y a , respectively, we can 
express the Yukawa couplings of down quarks and charged leptons at Mgut 



Y d = y GUT + k d j^Yv , Y t T = Y GVT + k e -£-y„ , (67) 

Mpi Mpi 

where a = 0(M G tjt)- The coefficients kd and k e are determined by the direction of the GUT breaking vevs. 
The relative transposition between Yd and Y\ is due to their embedding in SU(5) multiplets. Even though we can 
calculate the masses of the fermions at AfcuT with fairly good precision, we cannot fix the various couplings. The 
reason is simply that the observed mixing matrices diagonalise the products or combinations of Yukawa matrices. 
In the simplest case, where all matrices but Yd and Yi are diagonal, the quark mixing matrix diagonalises 

Y d Yj = roux^uT + k dJ ^ (FgutFJ + F CT F^ UT ) + Ud^j Y a Y^ , (68) 
while the leptonic mixing matrix diagonalises 



2 



YiY} = ^gut^gut + k e ^- C^ut^* + YjY£ VT ) + [ k ^J Y JY* . (69) 

(Again, these relations hold at Mgut-) In addition, as indicated by the factors k, the matrices are model- 
dependent. 

Since we do not want to restrict ourselves to a special version of a particular model, we proceed as follows. In 
the basis of diagonal charged lepton Yukawa coupling one gets 

V X = u\v d U 2 + j^-K (70) 



4 Here, we neglect the higher-dimensional operators which contribute equally to and Y;. In this discussion, we can absorb them 
in Ygut; however, they become important for B and L violating processes [66,67 . 
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SPS 


la 


lb 


2 


3 


4 


A 


B 




< 


0.0000032 


0.0000053 


0.000062 


0.0000045 


0.0000082 


0.0000103 


0.0000103 


\Sl 3 R 


< 


0.012 


0.019 


0.232 


0.017 


0.028 


0.036 


0.039 


\5f R 


< 


0.009 


0.015 


0.182 


0.013 


0.022 


0.028 


0.030 



Table 8: Upper bounds on \S}? R \, \^ R \ and \Sj^ R \ for the mSUGRA scenarios from lj — > Iff. 



with unitary matrices Ui and U<i. The starting point for universal A terms (cf. Eq. (|58|0 is the renormalisable 
Yukawa coupling, 

A l = A d = A Y GUT = A (y, - kej^-Y^j (71) 

Now we know that the entries of Kgut are generally of the right order of magnitude. Since the contributions 
from Y a are suppressed by a factor Mqut/A/pi & n d bottom-tau unification works well, they do not change the 
third generation's entries significantly. Then we can approximate the effect of the higher-dimensional operators 
with an additional rotation in the 12-sector, parametrised by one single mixing angle 9, 

(72) 

This parametrisation is similar to the one of Ref. |73| . We restrict ourselves to the lepton sector and place our 
boundary conditions at the GUT scale rather than the Planck scale, since we aim at constraints on the GUT 
parameters derived from low-energy data on LFV. In Section 14.31 we saw that the the LR off-diagonal elements 
of the slepton mass matrix are negligibly small. These elements are expected to become sizable now, due to 
inclusion of the additional mixing, parametrised by 9. In order to be as model-independent as possible, we will 
continue to assume universality of the soft SUSY terms at Mqut- Then the mixing does not affect the derived 
results for 5 % l L because the LL elements are not sensitive to 9. In a given GUT model, it may be more natural to 
assume universality at Mpi, as is the case in Refs. |33U72| . Naturally, their results are model-specific. 




Numerical Results for the LR Sector 



In order to derive upper bounds on 



8 V 

°LR 



, we assume all other off-diagonal elements are zero. The results for the 



different scenarios are listed in Table [H We show the relation between the branching ratio BR(/j, — > ey) and 8j^ R 
in Figure [131 The main contribution comes from a bino exchange which is independent of tan/3, contrary to the 
LL elements. 

Assuming diagonal slepton mass matrix at the GUT scale, the generated 8 l [ R at the weak scale depend basically 
on Aq and 9. The bounds on 8]^ R and 8 2 ^ R are too loose and the generated off-diagonal elements stay far below 
them. Only the 12 element can reach the experimental sensitivity. As long as one chooses Aq = , Sj^ R is 
negligible small even for a large mixing angle 9. 

Let us now vary Aq. This variation slightly modifies the mass spectrum at the electroweak scale via the RGE 
but the upper bounds on 5}^ R do not change significantly. For instance, in a modified SPSla scenario with Aq 
varying from —200 to 0, the bound lies within (3.22 — 3.34) • 10~ 6 . The generated S^ R element, however, can 
quickly exceed the experimental bounds, even for small values of 9 (Fig. I14p . Then we can derive a relation 
between Aq and the maximal allowed value for 9\ Figure [15] shows the maximally allowed value for 9 as a function 
of Aq for the SPSla and lb scenarios. The additional rotation reflects the different flavour structure of the down 
and charged lepton Yukawa couplings (see Eq. (|7"T|) ). Given the relation ([70]) . we conclude that for sizable Aq, 
the higher-dimensional operators respect the flavour structure of the tree-level couplings. 

In contrast to the LL sector, the results for the PMNS and CKM cases do not differ much for the LR sector. In 
both cases, the generated Sj^ R are negligible for vanishing mixing, 9 = 0°. We can easily understand this behaviour 
as we read off from Eq. (fBTj) that S]^ R contains the mixing from both Y v and Y\ . Any mixing coming from 9 ^ 
contributes to Y; and dominates over the mixing in Y v . Hence, there is no significant difference between CKM 
and PMNS cases. Ref. |76| has arrived at a similar conclusion, stating that intrinsic, arbitrary flavour violations 
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Figure 13: BR(fi — > erf) x 10 11 as a function of 5]? R for the different mSUGRA scenarios: From top to bottom: 
SPSla: red, SPS3: orange, SPSlb: green, SPS4: yellow, A: b light blue, B: brown, SPS2: blue, experimental 
upper bound black dashed. 



8.x 10 



6. x 10 



2.0 



^5 

■^s 0.000015 




Figure 14: Sj^ R as a function of 8 and the experimental bounds (red). Left hand side: SPSla with Aq — —100 GeV 
(blue dotted). Right hand side: SPSlb with A = -100 GeV (blue dotted) and A = -10 GeV (green dashed). 



in the slepton mass parameters at the high scale, even if relatively small, can completely obscure the loop effects 
induced by the seesaw mechanism. 



5. Conclusions 

Apart from neutrino oscillations, lepton flavour violating (LFV) processes have not been observed up to now 
and the individual lepton numbers have succeeded as good quantum numbers in charged lepton decays. Weak- 
scale supersymmetry, however, generically introduces an additional source of flavour violation. Hence, these rare 
processes enable us to study the supersymmetry breaking sector. In this paper we have performed a comprehensive 
study of the quantities S^y parametrising the flavour structure of the leptonic soft supersymmetry-breaking terms 
in the MSSM. Novel features of our analysis are the consideration of mass and anomalous magnetic moment of the 
electron and the (finite) renormalisation of the PMNS matrix by supersymmetric loops with soft terms. Further, 
we include tan /3-enhanced two-loop corrections to the LFV decays lj — > ^7 in a diagrammatic approach. Unlike 
previous analyses our method a priori does not involve any expansion in v 2 /M| USY , which becomes questionable 
in the case of large slepton mixing. We have subsequently expanded the exact result in w 2 /M| USY and have 
checked the accuracy of the expanded results. Our analysis of the PMNS matrix and the radiative decays follows 
the line of Refs. [Tg|-j20|, which have addressed similar problems in the quark sector. We finally analyse the effect 
of dimension-5 Yukawa couplings in the context of SO(10) GUT scenarios. 

Studying the one-loop renormalisation of lepton masses and PMNS elements at large tan /3 we have found 
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Figure 15: Upper bound for 9 as a function of Ao for SPSla (down blue) and SPSlb (top red) with varying Aq 
(in GeV). The allowed (Aq,9) region lies below the curve respectively. 

potentially large finite loop contributions to the electron mass m e while corrections to the PMNS matrix stay 
rather small. Applying a standard naturalness criterion to m e leads to the requirement that the loop contributions 
must not exceed the measured value. As a result we find |<5i£<5}fj?| ~ 0. 1 which involves the otherwise poorly 
constrained quantity S]^ R . The same parameter combination is constrained by the anomalous magnetic moment 
of the electron, a e , for which the MSSM contributions decouple if the corresponding mass scale Msusy becomes 
heavy. a e gives the same constraint on |^i<5)^| as m e for Msusy = 500 GeV. Further we have pointed out 
that the flavour-changing counterterms renormalising the PMNS elements generically appear in the charged-Higgs 
couplings, even if the latter are expressed in terms of weak neutrino eigenstates. The corresponding loop-corrected 
vertices are summarised in Eqs. (|3T j) (f33 | . Our two-loop corrections to the LFV radiative decays change the decay 
rates by up to 10% for large values of tan /3. Results on the inferred bounds on \S^ L \ for selected MSSM parameter 
points can be found in Table Q] Assuming reasonable SUSY masses < 1 TeV we find that BR(t — > ey) severely 
limits the size of the loop correction 5U e z to the PMNS element U e z\ For Msusy < 500 GeV we find \5U^\ < 10~ 3 
corresponding to a correction to the mixing angle #13 of at most 0.06°. Therefore SUSY loop corrections cannot 
fake a deviation from U e 3 = implied by tri-bimaximal neutrino mixing, if this PMNS element is probed with 
the precision of the DOUBLE CHOOZ experiment. Stated differently DOUBLE CHOOZ will probe the Yukawa 
sector and not the soft SUSY-breaking sector. 

The bounds on S^y are known to be severe, motivating the assumption that the SUSY breaking terms respect 
the SM flavour structure. As the symmetries and the particle content of the standard model point towards 
grand unification, one frequently assumes that these terms are universal at the scale A/gut, where the SM gauge 
couplings converge. Then the RGE running generates non-vanishing 8^ Y at the weak scale. In Sec. 4 of this 
paper, we have considered the Yukawa structure of two simple GUT scenarios. We calculated the size of the 
generated for various SUSY spectra, using the RGE for the MSSM extended with singlet neutrinos. The 
comparison with our bounds obtained from lj — > ^7 allows to constrain or even to exclude particular scenarios. 

In our RGE analysis we include the effect of higher-dimensional Yukawa operators (of dimension 5 or higher) 
which are needed to reconcile Yukawa unification with the experimental values of the fermion masses of the 
first two generations. If SUSY-breaking occurs above the GUT scale, flavour universality will naturally align 
the trilinear breaking terms with the dimension-4 Yukawa couplings, leaving the higher-dimensional terms as 
potential new sources of flavour violation. We have parametrised this effect by a new mixing angle 9 in Eq. (|72jl . 
For typical values of the universal trilinear term Aq one finds very stringent bounds on 9, as depicted in Fig. 1151 
As a consequence, the flavour structure of down-quark and charged-lepton Yukawa couplings must be similar 
for sizable Aq. This result hints at flavour symmetries which are respected by the higher-dimensional Yukawa 
operators. Note that it also applies to renormalisable couplings with a higher-dimensional Higgs representation, 
which couple differently to down quarks and charged fermions. In addition, the higher-dimensional Yukawa 
operators are generally consistent with all symmetries, hence appear naturally and yield significant corrections 
to the light generations' masses. The same qualitative result, aligned flavour structures of dimension-4 and 
higher-dimensional Yukawa couplings, has been found in a complementary analysis of the quark sector |75| . 

With the upcoming LHC experiments we will explore whether weak-scale supersymmetry is realised in nature. 
In addition, new flavour experiments like MEG will probe lepton number violation. Our analysis stresses once 
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more the importance of lepton flavour physics to map out the parameter space of the MSSM. Our GUT analysis 
exemplifies the well-known potential of lepton flavour physics to probe theories valid at very high energies. 
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A. Loop integrals 

We list the loop integrals, used in Section [3] 

B {x,y) = -A — hw —\n\ with A = - - lE + hi47r (73) 

x — y \ir y — x fi 4 e 

xy In - + x z In - + yz In ^ 

= (s -„)(*-,)(„-,) (74) 
, , \ fi(x,y,z) ~ fi{x,y,w) 

J2{x,y,z,w) = (75) 

z — w 

w In — y In ^ z In - 

Z i Z(j i X 



(w - x)(w - y)(w - z) (y - w)(y - x)(y - z) (z - w)(z - x){z - y) 
P , \ f2(x,y,z,v) - fo{x,y,z,w) 

F Q (x,y,z,v,w) = (76) 

v — w 

xy In - + xz In - + yz In - 

C (x,y,z) = - % y *, ' (77) 



(x - y)(x - z)(y - z) 



x y 2 In ^ z 2 In - 

ft^.-A-ta^-^JL^-^-^, ,78) 

y 2 ln^ z 2 ln^ w 2 In & 

D 2 (x,y,z,w) = — ^ \- --- \. r- ? ^ r, (79) 

{y - x){y - z){y - w) (z - x)(z - y)(z - w) [w - x){w - y){w - z) 

F?(x) = - 2 ,, \l-6x+ 3x 2 - 6x 2 logx] , 
(1 — xy 

F2{X)= (1^)3 t 1 - a;2 + 2a;l0 g a; ]' 

FP(x) = t- 2 [2 + 3x - 6x 2 + a; 3 + 6a; log a:] , 

(f — a;) 4 

^ (a;) = 2(1 -xf + AX ~ ^ ~ 2 bg ^ ' (80) 



B. Interaction of gauginos, sfermions and fermions 

The convention and notation of [38] is used with some little modification. The factors y/2 associated with the 
vacuum expectation value is omitted, such that 

v = \l v l + v l = 174 GeV ( 81 ) 

and the ratio between the vacuum expectation values are denoted by tan/3 = v u /v d . The Yukawa couplings are 
defined as follows 

m u = v u Y u , m d = -v d Y d , mi = -v d Yi ■ (82) 
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Neutralinos x.i 



* = (B,W,H° d ,H°) , 



L x° m , 



M N = 



--(*°) 
2 V ' 

( Mi 



31 "d 
' V2 

\ %/2 





92V d 

V2 



h.c. 




31 fd 


Sl"„ \ 




%/2 \ 


9'2'',/ 


32 fu 









-/i 


-A* 


/ 



(83) 



Mjv can be diagonalised with an unitary transformation such that the eigenvalues are real and positive. 

so 



Z n MnZn 



M 



(84) 







(M^) 2 can be used. Zn consists of the eigenvectors of the hermitian matrix 



For that purpose, 

M^ n Mn. Then the columns can be multiplied with phases e 8 ^, such that ZJ^M^Zn = Mjy has positive and real 
diagonal elements. 



Charginos xf 



W + ,HJ,W-,HJ , 



-i (4- ± ) T M c * ± +h.c. 



X 
X 



M 2 g 2 Wi 
52 A* 



(85) 



The rotation matrices for the positive and negative charged fermions differ, such that 





Z T XZ+ 



•"XI 

m 



\2 



(86) 



Sleptons and sneutrinos 

There are two ways of arranging the sleptons in a vector, either by family or by chiralities. The latter approach 
is adapted for the most general case and is used in Ref. [38] , whereas the former is convenient if LFC is assumed 
or for small off-diagonal elements treated as a perturbation. In this latter case, we have 



L, m - 2 



Mf 



A l R ~A P-R f l f R ) M f (CL e R Jl L m TL fR.) T 



( < 


m l LR 




Am T R 


Am LL 


A mT R \ 


9 

mi 


m l R 


Am RL 


Am RR 


Hi 




A <1 


Am i e R 




m kn 


Aml T L 


Am LR 


Am RL 


Am RR 


mf, 

URL 


9 

mi 

Mi? 


Hl 


Am RR 


Am T L e L 


A ™ T L e R 




Am LR 


m 2 TL 


m L 




A ™RR 


Am RL 


Am RR 


m r RL 


m l L R / 



(87a) 



with 



ml 


=<*> 






l) LL 





(-1) 



LJRR 



LR 



RL : 



4 s 2 r 2 



$ij + vjY l 2 S lj + (ml)ji 



2c 2 



(87b) 

(87c) 
(87d) 
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The rotation matrix Z is denned as 

Z^MfZ = diag (m?,...,mg). 
This matrix can be spit up into three parts 

(Z e \ 

(89) 

In general, these Zi are 2x6 matrices, which reduce to 2 x 2 matrices, respectively, with zeros in the remaining 
entries in case of vanishing LFV. Then for every generation one can write: 

°*V » = W (90) 



Alternatively, in the former case, one defines L\ := f£ and i? 7 := ej^, which mix to six charged mass eigenstates 
Li i — 1 ... 6, 

L 2 = Z L L i , R =Z L LJ, ZH 2 t 2 IZl = drag (m Li , ... ,m L J . (91) 

\\" l L)RL \' n L)RRJ 

These two approaches can be translated into each other by the following substitutions (with i = 1,2 for LFC 
and i = 1....6 for LFV): 

^li = Z H = = Z U ^4, = Z 2i Z 5^ = ^2, = ^2i _ /g 2 ) 

The sneutrinos and the left-handed sleptons have a common SUSY breaking soft mass. The weak eigenstates 
L\ = vi can be rotated in the sneutrino mass eigenstates i/j via Z v , 

L\ = Z»v s , ZlM\Z v = diag (ml ,m? 2 ,m? 3 ) , .M* = f!i?fz_?v) i + m | . (93) 

w c w 

Lepton-slepton-neutralino 

• incoming lepton I, outgoing neutralino and slepton If. 

iTf U = z (g lZ % + g 2 zi) + YtZpZ 3 /) P L + i (- 9l V2Z?z)i* + Y l Z\ i Z^*) P R 



• outgoing lepton I, incoming neutralino and slepton If. 

i (vf'J = i (- 9l V2Z, ! "Z« + Y,Zt«Z%) P L + i ( gi Zl" + S2 zy) 

Lepton-sneutrino-chargino 

• incoming lepton, outgoing sneutrino and chargino: 

ilf** = -i {g 2 Zl k P L + YuZ^Pr) Zr 

• outgoing lepton, incoming sneutrino and chargino: 

i (lT X ) = -* (Y h Z 2k P L +g 2 Z\ k *P B ) tfj 



R 
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C. Renormalisation group equations 

In the following fi denotes the energy scale (and not the 11 parameter of the superpotential) and t = ln(/x). For 
cji the GUT normalisation is used {gi UT = J\gf M )- 



At one loop order the gauge coupling in the MSSM evolve according to 

d . . 1 66 2 , . 
_ ai(t ) = __ ai ( t)> 

|a*(t) = --Uaft). 



(94) 



For the running of the gaugino masses, we use that gf(i)/Mi(t) is independent of the scale t at one loop order. 
Defining k = g\{tG\ST) /"ii/a an d assuming universal gaugino masses mu2 at the GUT scale, you can use Mj(t) = 
gf{t)/k in the RGE. 

The running of the Yukawa couplings and the Majorana mass matrix between Mqut and Mr at one loop level 
is given by the following differential equations: 



Ay 

alt u 
d 

dt Yd 



1 



d_ 

dt 
d_ 



5- 



16tt 2 
1 

16tt 2 
1 

16tt 2 
1 

16tt 2 
1 

8^ 



tr (3y„yj + F „ F t) _ K> fl a _ 3 ff 2 _ H ff 2 ) ! + 3^ + y Jy, 



1G 



7 



tr [3Y d Yj + Y t Y?) - fgi - 3 ff 2 - - ff 2 J 1 + 3YjY d + YjY u 
tr (3F„Fj + U f) _ 3ff 2 _ 3^ ! + 3Y }Y V + Y?Y t 

tr (sYdYl+Ytf 



3.92 - pi) 1 + 31?*j + 



M fl (y yj) + y u y}m r 



The running of the Aterms is given by 

1 



dt l 



16tt 2 



a u ^r(3yyj + yyj) - y ff3 2 - 3. 9 2 - 

( 32 26 

+ y ( tr(6yj A u + 2Y}A V ) + —g 2 3 M 3 + 6g 2 2 M 2 + -3^ 



dt 



-A, 



16tt 



+4yyjA u + 5A U Y^Y U + 2Y u YjA d + A u Y\Y d 

1 [ A d (tr(3yyj + yy/) - ^ - 3 ff2 2 - ^ ff 2 ) 

(32 14 \ 

tr(6y d U d + ZY^Ai) + —g 2 3 M 3 + Qg 2 2 M 2 + ^<7 2 Mi) 

+4yyJ A d + 5A d Y^Y d + 2Y d YjA u + A d Y*Y u 



dt 



-A. 



16tt 2 



tr(3y 1 yj + yyj)-3 ff2 2 -^ 2 



d A 

Jt Al 



16tt 2 



- y [tv{6Y^A u + 2Y}A V ) + 6g 2 M 2 + -g 2 M x 
-4Y V Y]A V + 5A v YjY v + 2Y V Y? A t + A v Y^Yi 



Ai ( tt(3Y d Yj + yy/) - 3. 9 2 - -g 2 

1 8 

+ y ( tr(6y]A d + 2Y?A t ) + 6g 2 M 2 + —g 2 ^ 



(95a) 
(95b) 
(95c) 
(95d) 
(95e) 



(96a) 



(96b) 



(96c) 
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■HYiY^At + SAitfY + 2Y l Y*A v + A t Y^Y v 



Sfermion- and Higgs masses (Notation: = nig, ml = to 2 ,, m 2 = m 2 , m| = m\, m| = m 2 ): 



with 



dt 



dt 



dt 



dt 



dt 



dt 



dt 



1 



Jt Q ~ Tte* 



IqYZY u + Y^Y u rriQ + m 2 Q Y\Y d + Y\Y d m 2 Q 

- 2 (Y d W d Y d + m 2 Hd YjY d + A\A d ) + 2 (Y^m 2 u Y u + m 2 ^!^ + At A u ) 



9 39 1 

15 ,2|)i,f |2 « „2 1 ?! r 1 2 oz „2|^ |2 i 1 



5f5 1 



^ [2 {mlY u Yl + YuYjml) + 4 {Y u m 2 Q Y^ + mf^Yj + 
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16tt 2 



16tt 2 



15 

2 (m^Yj + Y d Yjm£) + 4 (y^y] + m 2 Hi Y d Yj + A d A 



-Y^\ M ^-^9l\M 3 \ 2 + \g 2 S 



to|Y/Y + Y'Yml + m 2 L Y^Y v + YjY 



16tH 



1 



16tt 2 
1 

16tt 2 



16tt 2 



+ 2 (Y e ^m 2 e Y e + m 2 H YjY e + A\A e ) + 2 (Y^m 2 v Y v + m 1 H Y^Y l , + A\A V ) 
-^- g \\Mr\ 2 + %g 2 \M^-\g 2 S\\ 

2 (mlYiY? + YY/to 2 ) + 4 (Y e m 2 L Y^ + m 2 H Y e Y^ + A e A\) 
-| 5 2 |A/i| 2 + ^) 1 

[2 {mlY v Y} + Y„Yjml) + 4 (Y v m 2 L Yj + m 2 H Y„Y} + A„4)] , 
[6tr (Y^m 2 Q + ml + m 2 H t)Y u + A^A U ) 

+2tr (Yj(ml + ml + m 2 H \)Y„ + A\A V ) - ^M^ 2 - 6g 2 |M 2 | 5 

6 tr M (m 2 Q + ml + m 2 H t)Y d + A\ A d 



5^ 



+2tr y/(m| + m 2 e + m 2 Hd l)Y, + A\A X - -g\\M x \ 2 - 6g 2 \M 2 \ 2 - -g 2 S 



S = tr i 



™d ~ 2m u -m L + m e ) - m H + m 



(96d) 



(97a) 



(97b) 



(97c) 



(97d) 

(97e) 
(97f) 
(97g) 

(97h) 



(97i) 



The neutrino Yukawa coupling Y v decouples from the RGE below the Majorana mass scale and thus disappears 
from the equations. Some peculiarities occur if you integrate out the right handed neutrinos separately, as we do. 
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